mathematical modeling on the karkheh reservoir stresses and its application to the dalpari fault

in this research, a new mathematical modeling on strength changes due to reservoir elastic stresses across the preexisting fault plane is introduced. the method has been applied to the dalpari fault, which is one of the potential seismic sources in the vicinity of the karkheh reservoir. in this method the distribution of total stress across the fault cannot be determined because the initial stress is unknown; the pore pressure due to the reservoir is also not considered. the mathematical modeling method has been explained briefly in the following. the lake first is divided into small rectangles of sides a and b by two sets of orthogonal straight lines, one set conveniently east-west and the other north-south. the mean water depth h in each rectangle with area s is estimated, and the water pressure on the floor of the rectangle is replaced by a vertical force f =? gs h at the center of rectangle. it is clear that rather smaller rectangles lead to more precise modeling, hence, each rectangle with increasing h is divided into some parts. the water pressure of the lake is simulated by a set of point forces f which applied in the -x3 direction and acting on the rectangles. we define now a mathematical model of the single force f in the elastostatic fields using the delta function conception: the point force f is defined as: the component of displacement at point due to in the direction, , is given by: where j=3, for the water pressure of the lake, and is the distance from origin to point p. the general three dimensional relationships between nine cartesian strain component and three cartesian displacements are given by: these nine terms constitute the infinitesimal strain tensor, a symmetric tensor with six independent quantities. the stress tensor is given by stress-strain relationships based on constitutive law called hooke`s law is given by: using the conception of the stress tensor and well-known relationships in elastostatic theory the various stress parameters such as shear and normal stresses due to reservoir can be determined at the point p in a plane with normal n. in this way, we would be able to achieve the strength values due to the reservoir across the specified preexisting fault plane. the shear stress ( ) and strength ( ) due to the reservoir across the preexisting fault plane are, respectively, as follows. is angle measured in the preexisting fault plane between resolved shear stresses due to reservoir and ambient causes, and it is measured from the direction of the latter coefficient of friction along the fault plane, is coefficient of friction across the preexisting fault plane. the earthquakes near new reservoirs is classified into the following three cases on the basis of positive, negative and zero values of ; a reservoir based on this classification may stabilize some parts and destabilize other parts of the same nearby fault surface. case i: induced or reservoir assisted natural earthquakes; > 0. this situation will arise when 0 ? ? < 90 and the reservoir stresses have a net destabilizing influence on some parts of fault plane in which has a suitably large component in the same direction as initial tectonic shear stress, therefore the earthquake occurs earlier than its natural time. case ii: natural tectonic earthquakes despite the inhibiting influence of the new reservoir; < 0. this situation will arise when 90